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ABS TRACT 


The heat transfer characteristics of yawed cylinders was measured 
in an open induction tunnel in the subsonic range. Two model diameters 
were tested, 0.50 inches and 0.25 inches. The espa lle number based on 
cylinder diameter was varied from 2750 to 33200 for the 0.50 inch dia- 
meter model and from 1430 to 16700 for the 0.25 inch diameter model. 

The tests showed a good agreement for both models with the experi- 
mental results of other investigators for the normal case, or zero yaw 
case. As yaw angle was increased from the normal position to about 35 
degrees, a peaking in the heat transfer was obtained. As yaw angle was 
increased further the average Nusselt number decreased. This behavior 
is attributed to the end effect of the models used whose aspect ratio 
were limited by the size of the test section. At larger yaw angles 
there is an agreement of the heat transfer data for models with similar 


aspect ratio. 
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T. Introduction. 


The problem of the forced convection heat transfer in flows normal 
to cylinders has been investigated since the early workers in the field 
of the hot wire anemometry. The mathematical description of the rate of 
heat transfer as a function of the air velocity was first presented in 
1914 by Professor L. V. King fiji. The relation, King's Law, has the 
form 


=A. +B nea; 2 
Nu = 1 1 


The form of this relation was subsequently found to be correct for wide 
ranges of Reynolds numbers and is particularly accurate in the velocity 
range from 0 to 300 feet per second for ordinary hot wires. The heat 
transfer from a hot wire thus varies with the square root of the wire 
Reynolds number if the temperature and composition of the fluid remain 
constant. King also measured a reduction of heat transfer by yawing the 
cylinder axis from normal to the flow direction. 

McAdams [2] presents the data of Hilpert, which includes cylinders 
with diameters of 0.0079 to 5.9 inches normal to the flow direction, 
together with a recommended curve relating the Nusselt number with the 
Reynolds number. The entire range of data is segmented into portions 
correlated by equations in the form Nu = B Re". Table 10-3 of refer- 
ence [2] gives the value of the coefficient B and the exponent n for 
several ranges of Reynolds numbers. Hilpert's data are the backbone 


of the present normal correlation. 


1 e e e 2 
Numbers in brackets refer to numbered items listed in references. 
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Squire, cited by Schmidt and Wenner [3], presents a theoretical 
solution at the forward stagnation point of a cylinder. The data 
obtained by Schmidt and Wenner verified precisely Squire's correla- 
tion as Nu = 1.01 eee 

Baldwin, Sandborn and Laurence [4], presented a summary of the data 
for the normal case in the continuum, slip, and free moleculae air flows. 
Giedt [5] investigated the effect of turbulence of the incident air 
stream on local heat transfer on a cylinder. 

The effect of sweep-back on the heat transfer has been largely in- 
vestigated, because of its application in the flight of aircraft and 
missiles, at supersonic speeds. Blunting the leading edges of wings 
reduces the local increase in heat transfer, but increases the drag. 
However, if the leading edge is swept back, the heat transfer rate due 
to the leading edge is reduced below that occurring if the leading edge 
is normal to the direction of flight and also the drag due to the 
blunted leading edge is reduced. References [6] and [7] investigated 
the effect of yaw on the heat transfer coefficient at supersonic speeds. 

The yawed infinite cylinder simulates approximately the leading 
edge of a swept-back wing and allows a basic simplification of the 
boundary layer theory. Sears [8] and other investigators observed that 
for incompressible flow over a yawed infinite cylinder the boundary layer 
development in the chordwise direction (normal to the cylinder axis) 
is independent of the spanwise flow. For compressible flow, however, 
this independence principle does not apply. 

Where the independence principle applies, the solutions for the 
boundary layer development in the chordwise plane are those which have 


been obtained for incompressible two dimensional flow. Sears presented 
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a series type solution to the spanwise velocity gradient distribution 
around a cylinder. It was noted [8] that in a steady, incompressible 
flow the distribution of spanwise velocity is the same as the distri- 
bution of temperature when the cylinder surface is maintained at a uni- 
form temperature different from the stream temperature, provided that 
the Prandtl number is unity and viscous heating is neglected. Based on 
this, Goland [9] showed that in an incompressible fluid of Prandtl 
number unity the heat transfer of a yawed cylinder varies as the square 
root of the Reynolds number based on the normal component of the stream 
velocity. Thus, for a given stream velocity, the heat transfer co- 
efficient decreases as the square root of the cosine of the yaw angle. 
This decrease is associated with the increase in boundary layer thickness 
due to yaw. 

The fact that yawing of cylinder reduces the average heat transfer 
coefficient has been shown by workers in hot wire annemometry. Reshotko 
[10] and Reshotko and Beckwith [11] gave a solution of the compressible 
laminar boundary layer equations at the stagnation line of yaw cylinders. 
For Mach numbers from 0 to 0.5 they found that the effect of yaw can be 
accounted for by writing Nusselt number as a function of the square root 
of the product of the Reynolds number and the cosine of the yaw angle. 
Thus, for incompressible flow the theoretical variation of the ratio of 
the stagnation line heat transfer for yawed infinite cylinders to the 


Stagnation line heat transfer coefficient at zero yaw reduces to 





h 
= - Guern 
/\=0 


where /\ is the angle of yaw, i.e., the angle between the spanwise axis 


of the model and the plane perpendicular to the free stream flow direction. 
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Thus a cylinder normal to the flow is at zero yaw, while the yaw angle 
for a cylinder parallel to the flow is 90 degrees. This was obtained 
experimentally with hot wires by Schubauer and Klebanof [12]. The 
Reshotko and Beckwith laminar theory cannot be applied to large yaw 
angles, since the initial hypothesis in the development of the theory 
that the spanwise derivatives are identically zero does not apply. 
Based on Reshotoko and Beckwith investigation, Baldwin, Sandborn and 
Laurence [4] indicated that for yaw angles between 0 to 70 degrees, 
the use of the correlation of convective heat transfer for transverse 
eyiinders, tit they presented in their Figure 6, with Re Cos/\ in- 
stead of Re, should give a good engineering estimate of the heat trans- 
fer for yawed cylinders for Re Cos/\ )400 in subsonic flow. 

Schubauer and Klebanoff [12] suggest that the yaw data could be 


correlated by the normal Reynolds number as 


Nu = A, (Re Cos A vile + B 


1 
Sandborn and Laurence [13] indicated that in general this correlation 
is adequate only near zero Mach number and angles of yaw less than 70 


degrees. They correlated their data for all yaw angles within ap- 


proximately ten percent with an equation of the form 


Nu = [A, (Re Cos asp win + B,] Cos/N\ +[A, (Re Skid + B,] SinN 


obtained by combining the heat loss of wires normal and parallel to the 
flow. However, they found systematic variations rather than random, 
indicating a need for a still more accurate equation. 

The object of the present investigation was thus to determine experi- 
mentally the heat transfer characteristics of yawed cylinders for yaw 


angles 0 to 90 degrees in the Reynolds number range 1500 to 33000. 


iS 


1s Description of the Equipment. 


a. Wind tunnel. 

The tests were performed in an open induction tunnel shown in 
Figure 1. The suction side of four feet square is covered with a fine 
screen to reduce large scale turbulence. The settling section is two 
feet long with a cross sectional area of one square foot. 

A nozzle with a contraction ratio 6:1 connects the upstream 
side to the test section. An eight foot long duct connects the test 
section downstream to the suction side of a fan propelled by a 0.75 
HP AC electric motor. A gate valve at the exit of the fan controls the 
mass flow rate. 

b. Test section. 

The plastic test section is two feet long and has a rectangu- 
lar cross section of four by six inches. Figure 2 shows the test sec- 
tion. Two circular openings were made on the ceiling. The first one, 
located 8.5 inches from the upstream flange, was used to locate a 0.25 
inch Flow Corporation pitot-static tube to measure the free stream 
velocity. This vill be referred to as station No. 1 or the upstream 
Station. VTne ‘cond opening, seven inches downstream of the first, was 
used to locate the models tor the heat transfer tests. This will be 
referred to as station No. 2 or the downstream station. The center of 
this station is located 1.75 inches from the nearest lateral wall in 
order to perform the tests for all yaw angles from 0 to 90 degrees, 
that is, changing the position of the models from a normal position with 
respect to the flow, to the parallel position. 

A plastic circular cover with a small hole for the location of the 


pitot tube was used at station No. 1. A similar cover was used at 
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Station No. 2 to be placed instead of the model assembly during the time 
when the model was heated outside the test section. Two quick action 
spring locks press the covers against the test seeettaiese ceiling to mini- 
mize leakage. Static pressure taps were placed in the side and bottom 
walls. 

c: Models. 

A protractor from a drafting machine mounted in a six inch 
circular plastic piece, was used to locate the models at the desired 
yaw angle in a fast way and with an accuracy estimated of the order of 
one degree. 

Two models were used for the present tests. The first is a 
0.50 inch diameter cylinder shown in Figure 3. A 0.50 inch diameter 
stainless steel holder piece was mounted in the protractor. The model 
was made using three similar copper cylinders 0.75 inches long and 
separated by small bakelite insulating cylinders 0.1875 inches long. 

A special technique was used to get a uniform model without 
discontinuities and good finished surfaces. The bakelite cylinders 
were made of larger diameter and with the ends like rings to hold the 
copper cylinders in the correct position allowing a good alignment. 

The cylinders were bonded together using epoxy adhesive. After dry- 
ing, the assembly of cylinders was machined down to 0.50 inches diameter. 

The copper cylinder near the holder piece and the bakelite 
insulators located at its sides have a 0.062 inch diameter hole in the 
center to allow the passage of a 30-gage copper-constantan thermocouple 
wire to the boundary of the center solid copper cylinder or the model 
proper. The basic idea of this design of the model, using the two 
copper cylinders similar to the model proper at its sides is to mini- 
mize the effect of conduction in the axial direction. Figure 6 shows 


a sketch of the 0.50 inch model. 
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The length of the model was governed by the width of the test 
section. As the yaw angle was increased, longer models were possible. 
This was accomplished by replaceable bakelite sections which were at- 
tached. 

Three removable extension pieces were used. One 0.656 inches 
long was used in the 30 and 40 degree yaw position, one 1.469 inches 
long in the 50 degree position and one 3.969 inches long in the 60 to 
90 degree positions. 

The second model shown in Figure 4, is similar to the first 
but 0.25 inches in diameter, all the other dimensions being the same. 

d. Heater. 

The heater was made by coiling a 18-gage, nichrome wire covered 
with ceramic bead insulators around a steel cylinder. A second steel 
cylinder was placed concentrically around the heater coil and covered 
with asbestos tape. Both heater ends were covered with insulation to 
decrease the heat losses. Openings of 1.25 inches diameter were made, 
to allow the passage of the models to a position in which the model 
proper was in the center of the heater. A wooden frame supports the 
protractor mount during the heating. 

The power was supplied by a Powerstat variable transformer. 
The current and voltage required for the heater were measured by Weston 
AC meters. 

e. Instrumentation. 


Figure 5 shows the general disposition of the instrumentation. 


rod 


ig Manometers. 

A manometer board with a set of Ellison inclined 
and vertical manometers was used for the pressure measurements. The 
accuracy of these manometers is within one half percent of the full 
scale. 

The following manometers were used: 

(1) Ellison inclined manometer, 0-0.2", smallest 
division 0.001". 

(2) Ellison inclined manometer, 0-2", smallest 
division 0.005". 

(3) Ellison inclined manometer, 0-4", smallest 
division 0.010". 

(4) Ellison inclined manometer, 0-12", smallest 
Aggpieas.ong 0.020". 

(5) Ellison vertical manometer, 0-17", smallest 
division 0.100". 

(6) Ellison vertical manometer, 0-36", smallest 
division 0.100". 

Zs Thermocouples. 

Thermocouples were made from 30-gage copper-con- 
stantan wire. The junctions for the model and for the free stream 
temperatures were first calibrated with a Leeds and Northrup precision 
potentiometer with ice and steam as the reference temperature. The free 
stream thermocouple junctions were located 22 inches upstream of the test 
section. 

ae Recorder. 


The temperature difference between the model and the 


Ze 


free stream was registered during the experiments using a 7100B Hewlett 
Packard strip chart recorder with plug-in 17500/A. 

The free stream temperature was measured for each 
test with a Leeds and Northrup portable potentiometer with an ice refer- 


ence. 
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S Experimental Procedure. 


a. Test section calibration. 

The test section was calibrated without the models. The velo- 
city profile was determined at the two stations of the test section by 
a traverse with the pitot tube in both vertical ead horizontal direc- 
tions. At the same time static pressures were measured at station No. l 
and at the pressure taps corresponding to the position of the model at 
station No. 2. 

The velocity was found to be practically uniform across the 
tunnel test section to within 0.25 inches of the wall. Comparison of 
the velocity measured with the pitot tube and that evaluated by static 
wall pressure measurements, indicated that the velocity could be evalu- 
ated using the pitot static tube with a stagnation point pressure co- 
efficient of 1.0. Velocities from 11 to 145 feet per second were mea- 
sured at station No. 2. 

b. Heat transfer experiments. 

The models were carefully polished before each run to minimize 
radiation losses. The Powerstat was set for an output to the heater of 
11 volts. The heater was maintained operating at a constant current of 
4.5 amperes. 

The millivolt potentiometer was standardized before each run. 
The strip chart recorder was calibrated by supplying a voltage of one 
millivolt from the potentiometer. This was done first at the one milli- 
volt scale of the recorder and then at the two millivolt scale. 

Before each run the atmospheric pressure indicated by a Princo 


Fortin Type barometer was observed. The free stream temperature was 
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measured for each run by a thermocouple and the portable potentiometer 
and also by a mercury thermometer located at the inlet of the tunnel. 

The tunnel was set at the desired velocity by adjusting the 
control valve. Once the conditions of steady flow were met, the model 
was removed from the test section and placed into the heater. A plastic 
cover was used to replace the protractor-model assembly. 

The models were heated to about 65 to 105 degrees above ambient 

temperature, corresponding to a difference of 1.5 to 2.5 millivolts of 
thermocouple output. It took about twenty minutes for the heater interior 
steel surface to reach a temperature of about 250 degrees Fahrenheit. 
The heating process was accomplished by free convection and radiation. 
It took about five minutes to heat the 0.50 inch model to 140 degrees 
Fahrenheit and about three minutes to heat the 0.25 inch model to 150 
degrees. 

Once the desired model temperature was reached, the model was 
removed from the heater and positioned in the test section. The plastic 
cover piece of the model assembly was quickly fastened by the two lateral 
spring locks. The model was placed in the desired yaw angle using the 
protractor mechanism. 

The temperature of the model during the initial phase of the 
cooling process was monitored by the recorder on the five or two milli- 
volt scales as convenient. Just before a temperature difference corres- 
ponding to 1.0 mv. was reached, the recorder scale was set at this 
value and the temperature-time history of the model initiated at a con- 
veniently selected chart speed. Figure 7 shows a typical model cooling 
record. 

The manometers were observed during the cooling of the model and 


the readings noted when the temperature difference between the model and 
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the free stream indicated by the recorder was about 0.6 millivolts. In 
general the manometer fluctuations during the tests were small. Negli- 
gible fluctuations were observed at the low and medium Reynolds number 

runs for both models. At the higher velocity runs, fluctuations of the 
order of 1.5 to 2.0 percent of the reading were observed. 

Ten series of runs at different flow rates were made with the 
0.50 inch model. Seventeen runs were made at each flow rate with this 
model. The first ten runs were made using the basic model without ad- 
ditional cylinder extensions at the end. These runs were performed for 
each ten degrees of yaw angle from 0 to 90 degrees. Another seven runs 
were made using the additional cylinder extensions. 

The procedure followed in the individual runs with the 0.25 
inch model was the same except that only nine mass flow rates were 
used. Some vibration was detected during the two higher speed runs 
with this model using the largest extension. 

The pitot static tube was located 0.5 inch from the upper wall 
of the test section to avoid disturbances in the flow upstream of the 
model. For the runs with the 0.50 inch model the upstream velocity 
changed from 11 to 136 feet per second. The corresponding changes in 
Reynolds number was from 2750 to 33200. For the 0.25 inch model the 
velocity changed from 12 to 138 feet per second. The corresponding 
change in Reynolds number was from 1430 to 16700. 

It was noted that by changing the yaw angle of the 0.50 inch 
diameter model from the normal to parallel position, the air velocity 
increased about five percent at high flow rates, decreasing to about 
one percent at the low flow rates. However, for the 0.25 inch diameter 
model the increase in local velocity is less than two percent at high 


flow rates. 
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4. Presentation of Results. 


The experimental data is given in both tabular and graphical form. 
Tables III and IV give the heat transfer data for the 0.50 inch and 0.25 
inch models respectively. Figure 9 is a summary of the normal heat trans- 
fer data for both models. Hilpert's correlation, mentioned by McAdams 
[2], and data from Schmidt and Wenner [3] is also plotted for comparison 
purposes. 

Figure 10 shows the average Nusselt number versus Reynolds number 
correlated in two ways for the 0.50 inch model at 20 degrees yaw. 

Figures 11 through 13 show the average Nusselt number versus 
Reynolds number for both models at yaw angles 30, 50, and 60 degrees 
respectively. Figure 14 shows the parallel heat transfer data for both 
models. 

The experimental variation of the average Nusselt number for the 
yawed cylinder to the average Nusselt number of the cylinder at zero yaw 
at representative Reynolds numbers is given in Figures 15 through 17 for 


the 0.50 inch model and in Figures 18 through 20 for the 0.25 inch model. 
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oe Discussion of Results. 


a. Normal case. 

Figure 9 presents a comparison of the data obtained at zero 
yaw for both models. Hilpert's correlation cited by McAdams [2] and 
the data of Schmidt and Wenner [3] are also shown. 

In general the present experimental points lie above Hilpert's 
correlation, the maximum deviation being about 18 percent at the low 
Reynolds numbers decreasing to about two percent at the higher Reynolds 
numbers. The average of local measurements made on big cylinders by 
Schmidt and Wenner lie above the present experimental points and also 
above Hilpert's correlation. The maximum difference between Schmidt and 
Wenner experimental averages and the present tests is of the order of 
13 percent. Schmidt and Wenner explained the fact that their data lie 
above Hilpert's correlation to inevitable and greater measuring errors 
in the experimental method that they used. 

A correlation in the form suggested by Richardson [14] was 
obtained for the two models in normal flow. The mean heat transfer co- 
efficient for a cylinder can be estimated as the sum of the contributions 
from the forward boundary layer and from the rearward separated regions. 
At Reynolds numbers for which an attached turbulent boundary layer is 
absent, the mean transfer coefficient can be written as 


Nu = Cc, Rei /2 + C, re2/3 


a¥E » against Rel!® gives a 
Re” 
ie i. 
ed for the present normal tests for 1400¢ Re< 33200 is 
1/2 zy 





The graph of the mean heat transfer, as 


straight line of slope C, and intercept C The correlation thus obtain- 


Nu = 0.400 Re + 0.354 Re 


with a maximum deviation of eight percent. 
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The intercept obtained is in agreement with the expected value 
from the boundary layer theory. The slope is within four percent of the 
value given by Jakob and Hawkins [15] as the average from measurements 


of different investigators. 


b. Yawed cylinders. 

Reshotko and Beckwith [11] give a solution of the compressible 
laminar boundary layer equations at the stagnation line of yawed cylin- 
ders. For Prandtl number equal to one and the incompressible flow case 
their solution indicates that the ratio between the heat transfer co- 
efficient at an angle of yaw to the normal (zero yaw) heat transfer 
coefficient at the stagnation line is equal to the square root of the 
cosine of the yaw angle. This theoretical solution is also shown in 
Figures 15 through 20 for comparison with the present experimental 
data at different yaw angles, for both the basic model and for the model 


with extensions. 


The decrease in the average heat transfer coefficient for the models 
with extensions is apparent from the figures. The reason is due to the 
greater aspect ratio L_/D of the models with additional cylinders. 
These more nearly approach the ideal case of an infinite long cylinder, 
the basic assumption of the theory. The length of the models of the 
present tests was limited by the width of the tunnel test section. 

Table I presents the aspect ratios of the models for the different 
yaw angles. The aspect ratio is considerably smaller than the 0.5 inch 
diameter 24.0 inches long copper model used by Bechwith and Gallagher 
(Ah. 

The data obtained from the 0.25 inch model is closer to the laminar 
theory because the aspect ratio for this model is twice that of the 0.50 


inch diameter model. 
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The trend of the data is better shown in Figures 15 through 
20 which present the ratio of the average Nusselt number at a given 
yaw angle to the average Nusselt number at zero yaw, as a function of 
yaw angle for a particular approach Reynolds number. Also shown is the 
theoretical laminar solution from Reshotko and Beckwith [11]. 

For the 0.50 inch diameter basic model, increasing the yaw angle 
from 0 to 35 degrees, the average heat transfer coefficient increases by 
about 26 percent at a Reynolds number of 31900, see Figure 17. This trend 
is different from that indicated by the laminar theory and which was 
obtained experimentally for wires by Schubauer and Klebanoff [12]. 
Sandborn and Laurence [13] observed a heat transfer coefficient at 35 
degrees that was about ten percent lower than the values at zero yaw. 

Further increase in yaw angle from 35 to 70 degrees results 
in a reduction in the average heat transfer coefficient of the 0.50 inch 
basic model of about 24 percent from the peak values. The trend of the 
experimental data obtained at other Reynolds numbers with this model is 
similar, see Figures 15 and 16. In general at lower Reynolds numbers 
the peaking values are smaller. The reduction in average heat trans- 
fer coefficient for the tests with the extensions is indicated in these 
figures. 

The data obtained from the 0.25 inch model shows a similar 
trend with an increase in the average heat transfer coefficient of about 
12 percent at a yaw angle of 30 degrees at a Reynolds number of 9050, 
see Figure 20. An increase in the yaw angle to 70 degrees results in 
a decrease of about 38 percent from the peak value. The percent of in- 
crease in the average heat transfer coefficient at the peak value is 
smaller for the 0.25 inch model than for the 0.50 inch model. The maxi- 


mum increase in the average heat transfer coefficient for the 0.25 inch 
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model with the extensions is about the same percent as the basic model, 
see Figures 18 through 20. The peak in average heat transfer is at 
about the same yaw angle. At larger yaw angles the reduction in heat 
transfer coefficient of the models with the extensions is greater because 
of the larger aspect ratio. 

The trend for both models is similar to that obtained for 
Beckwith and Gallagher [7] in supersonic flow. They explained their re- 
sults by a phenomena of transition from laminar to turbulent boundary 
layer. This transition apparently occurred primarily as a result of a 
dynamic type of instability which is associated with the secondary flow 
in the boundary layer. The same agrument cannot be used to explain the 
trend of the present tests, because the range of Reynolds numbers is 
well below the values indicated by Owen and Randall, mentioned in [7], 
for thé transition from laminar to turbulent as a result of secondary 
flow that occurs within the boundary layer on a yawed cylinder. It is 
interesting to note that the Reynolds number for the transition phe- 
nomena described by Owen and Randall decreases with yaw angle. Burnsnall 
and Loftin, mentioned in [7], observed a decrease in the stream Reynolds 
number for transition on circular cylinders (Re,,) » from 3.7 x 10° to 
e.i” x 10° for an increase in yaw angle from 0 to 60 degrees. Fur- 
thermore, the experimental transition Reynolds number obtained by 
Scott-Wilson and Caps, mentioned in [7], was about one tenth of the value 
predicted by the criterion of Owens and Randall. Being the trend of 
the data obtained similar to that of Beckwith and Gallagher [7], how- 
ever, the peak values obtained in the present tests are well below those 
obtained by them. The maximum deviation from the value predicted by 
the laminar theory for the 0.50 inch model is about 44 percent, and 


is about 24 percent for the 0.25 inch model. The use of finite 
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length cylinder with an exposed end is the departure from the infinite 
length idealization of the theory. The pattern of the flow is changed 
from that corresponding to an infinite yaw cylinder, due to the end ef- 
fect of the model. The decrease in heat transfer coefficient obtained 
for the tests with the extensions was less than expected for the corres- 
ponding increase in the aspect ratio. The heat transfer characteristics 
of the models used in the present tests seems to be influenced by the end 
effect whose intensity depends upon the distance from the model proper 
to the leading edge and the yaw angle. This effect takes place as soon 
as the model is given a small yaw angle. For a given model at a given 
yaw angle the heat transfer coefficient decreases by increasing the 
distance from the leading edge to the model proper. For a model of a 
given aspect ratio the influence of the end in the flow past the cylin- 
ders increases with the yaw angle. At low yaw angles, the large aspect 
ratio of the 0.25 inch model resulted in values closer to the predicted 
than the lower aspect ratio of the 0.50 inch model. At large yaw angles 
the influence of the aspect ratio in the intensity of the end effect can 
be seen by comparing the heat transfer data of the 0.5 inch model with 
the larger extension with that of the 0.25 inch basic model both of 
which have about the same aspect ratio. The data for these two models 
is approximately the same. Change in aspect ratio results in a departure 
from this observed agreement. 

Cus Parallel case. 

The flow geometry for this case becomes completely different 
from the previous cases. The influence of the aspect ratio in the in- 
tensity of the end effect can be appreciated in this limiting situation. 
A big difference in the heat transfer data for models with different as- 


pect ratio can be noted. The agreement in the heat transfer data for 


312 


models with similar aspect ratio is seen in Figure 14. Thus data for 
the 0.50 inch model with the larger extension and the data for the 0.25 
inch basic model can be correlated in the form indicated by Sandborn 


and Laurence [13] by Nu = 0.48 =>. with a maximum deviation of about 


12 percent. 


$3 


6, Conclusions. 

The tests showed a good agreement with the experimental results of 
other investigators for the normal or zero yaw case. 

The peaking in the heat transfer coefficient obtained at a yaw 
angle of about 35 degrees is attributed to the end effect of the models 
used whose aspect ratios were limited by the size of the test section. 
This end effect causes a change in the flow pattern with respect to that 
of an infinite cylinder. 

For large yaw angles the agreement of the heat transfer data obtained 


for models with similar aspect ratio is apparent. 
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FIGURE 7. Typical Model Cooling Record 
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Average Nusselt number versus Reynolds number at 30 degrees 
for the 0.50 inch and the 0.25 inch models. 


Figure ll. 
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Average Nusselt number versus Reynolds number at 50 degrees 
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Figure 12. 
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Summary of parallel heat transfer data for the 0.50 inch 


Figure 14. 


and the 0.25 inch models. 


++ |} 
+ 
-fe 
=~ + 
Coe alld 


<< 


Serackacs: un al apie a8 HHH tH HEH 
Beane Shbes Cage sobs segesdebslugege H Bit BRE 
t ++ sae Huu Ge 4 t a asec. ag - 
Sb pe ii Lp SERS 
‘ sisi ae sags iat ae e AB SB ° 
; a ' eRe : iL Bis anne. Sa htt veeemncecgoeeee fe 
e 5 @tede r) sagees 

| ses tae of eserteze eae estat 

ty ° 7. eee 
TO 
2 pia i 5 4 28 at PEE 
| ee Cee Hd Utedes neoeem@ees ce ae 
a at yee: C4 sprite 


a a 
tees = 


C 
o sud inch model with extension <=: 


ce 


Hi 





, 
re 


He a 
sl a 


— . 


. ary . 
See tH 
f cil % 


eG Fo pg RR Ti ke aiess roth 
re soe ae ‘ THY eee mn 
Bs = ne 4 se t ais HE = 
: an <r oag 
g ae see ats Se ae Bete 
Sst > Egec! seeps 000g cu He HE a 


oor | 
ars a § 


EE: eee ea ee 
PCLCe eee 
EYE AT SHES Ye A ee SE 
PEP ea ages ie 
ee egores 44 seog +++ 

SID) Grae aaa 


emia ae 
Tore Fila. Pee 
See e e2e TEL 


wr 9” TER : poss ores: SESS! ISOS Sg 

* aed ao o he oom ee 
. ee o@ee ae no e- —_ 
4 bd . e e * eee eer - 

° . ee . 8. 2-7 = Sted 

ao ' @ree 

eee owe ee Cue es 

° eooeve . 

fe) oe ° ° e 

(P : 

. . ees ° - ° eee . rere iS 


etee 












e . e [os . 
ee é a we ° + eere 
: 
—. Co oe Ce -- o+ eeee eo b+: 
oe a 6 oe ot § eo oe a ee = ~_ eo 
> + Os dint Pm ed ote a leaps fatal seat 
ry ° ees ++ ee eee eee: —_ —, __ + + =c-) 
sd 
‘4 oe. - hag 1 -_ ooo swe eve . ° _ ° fatmta) 
° Cra} ' eee aoe oy ° . 0 $+ +e: 
J. -_ on far ae] 
eee ‘e ee ee = . 2 4 ewe 
ae — aa one allen oslo 
i MDG: tt? © eG eee aeae 
~e awert + os orbs oad pula eee 
elSesiwe hen = = Eeaa ee 
ee nae ¥ sas 
ae a . ta 
or ae sa razal 
an - (ate sees 
sept SE e 
erie = 74 aes 
Comes OG BGSQGR2S — es - ! meat al +4 toe 
gone g ee ens roses Das -T2}- “4 = eos ss a ese 
nooga fa SPegas resee pages sane) pie) eee OSs ee we $e mot se 
. eatin sian 
_ : <= i 





oad (ape coed 





“+? 





wy oe) = re a) aa wy 
° e e e e e 
= a =~ © © © 


Sl 


» degrees 


Yaw angle, A 


Reynolds 


2800 
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Figure 15. 
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» degrees 


Yaw angle,A 


Experimental variation of the average Nusselt number for the 0.50 inch yawed 


model to the average Nusselt number of the model at zero yaw. 


Figure 16. 
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Experimental variation of the average Nusselt number for the 0.50 inch yawed 


model to the average Nusselt number of the model at zero yaw. 


Figure 17. 
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Reynolds 


: 1460 


Experimental variation of the average Nusselt number for the 0.25 inch yawed 
Number 


model to the average Nusselt number of the model at zero yaw. 


Figure 18. 
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, degrees 


Yaw angle, /A\ 


Reynolds 


Experimental variation of the average Nusselt number for the 0.25 inch yawed 


model to the average Nusselt number of the model at zero yaw. 


Figure 19. 
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Table I. 


Diameter D in 
Length L in 
Mass m lbm 


Basic Length L in 


Aspect ratio L/D 


Data of Test Models 


Basic 


Models 


0.50 


Ow 75 


0.04755 


3.94 


7.88 


O225 
OA Tis, 
0.01188 
3.94 


15.76 


Aspect Ratio L/D of Models with the Extensions 


for the Yaw Angles Tested 


Yaw Angle 
Total Length 


Aspect ratio of the 
0.50 inch model 


Aspect ratio of the 
0.25 inch model 


30-4 


4.44 


8.88 


Li5 70 


0 


a) 


50 


De25 


10.50 


21.00 


60-90 


TRIAS, 


15.50 


31.00 


Table II. 


Material Density 
Ton pee 

Copper 25s; 

Bakelite 79.5 


Properties of Copper and Bakelite 


Specific 
Heat 


Btu/lbm°F 
0.0918 


0.40 


Thermal 
Conductivity 


Btu/hr £t°F 
225 


0.1088 


Thermal 
Diffusivity 


ct- fae 
4.352 


0.0034 


Properties of pure electrolitic copper were taken from Eckert 


[17]. The specific heat was estimated for t_ = 90°F. 


Properties of bakelite were taken from the Modern Plastics 


Encyclopedia [20]. 
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Table 


Reynolds Number 


Zh 
3654 
5264 
7461 
11487 
14137 
16531 
20671 
25932 
31380 


2749 
Bo ee, 
5265 
7462 
11420 
14042 
16448 
20463 
25486 
30909 


2802 
3543 
5290 
7463 
11418 
14043 
16363 
20530 
25596 
30909 


III. Summary of Heat Transfer Data for the 0.50 


Inch Diameter Model 


Basic Model 


Yaw Angle = 0 degrees 
(Normal case) 


Yaw Angle = 10 degrees 


Yaw Angle = 20 degrees 
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Average Nusselt Number 


S073 
2.0 
39.4 
50.4 
64.2 
70.6 
ho) 8. 
79.5 
93.3 
104.3 


30.8 
37.6 
40.6 
72 8 
64.5 
75 .0 
ie 
87.6 
97.6 
112.9 


Bo. / 
oy. 5 
45.1 
59.4 
75.1 
81.8 
91.0 
97.7 
ibe 
124.0 


Table III. 


Re 
2797 
3542 
5290 
7463 
11399 
14091 
16448 
20726 
25707 
31211 


2804 
3542 
5290 
7463 
11401 
14091 
16531 
20790 
25971 
31384 


2846 
3542 
5290 
7463 
11517 
14188 
16655 
20920 
26439 
31861 


2846 
3542 
5290 
7462 
11517 
14280 
16771 
21298 
26796 
32280 


(Cont) 


Yaw Angle = 30 degrees 


Basic Model _ 


Nu 
3453 
38.1 
46.7 
59.4 
76.0 
84.4 
94.7 
102.3 
122.8 
132,3 


33.6 
37.9 
46.7 
58.2 
75.8 
83.8 
Loco 
103.4 
12223 
131.3 


31.8 
Sais) 
44.6 
55.1 
70.6 
LOU 
87.3 
Gee] 
ey. 
124.9 


28.8 
3253 
3955 
50.6 
65.0 
72a) 
78.6 
S7-2 
101.1 
114.6 


Yaw Angle = 40 degrees 


Yaw Angle = 50 degrees 


Yaw Angle = 60 degrees 
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Model with extension 


Re 
2849 
3546 
5321 
7358 

11391 

14033 

16491 

20644 

25644 

31003 


2848 
3546 
5296 
7465 
11449 
14150 
16449 
20900 
25802 
31132 


2846 
3548 
5306 
7465 
11508 
14176 
16521 
20837 
26220 
31562 


2845 
3548 
5302 
7465 
11509 
14176 
16655 
21090 
26527 
32166 


Nu 

33.6 
35.6 
we I 
Sia 
70.3 
78.6 
88.2 
975 
115.6 
13228 


3330 
36.1 
45.0 
Byio 
72.8 
81.8 
88.6 
100.9 
Ll7a? 
125.7 


30.1 
34.4 
41.9 
op ee) 
66.4 
74.4 
83.4 
a2 
109.5 
Li9e2 


202 
285 
Sie 
45.7 
56.5 
62.9 
69.9 
77.3 
70.3 
102.8 


Table III. (Cont) 
Yaw Angle = /0 degrees 


Basic Model - Model with extension 
Re Nu Re Nu 
2846 Zon 2 2845 23.0 
3543 2325 3549 26. 7 
5315 34.0 5302 B57 
7461 43.0 7466 3738 
11572 55.6 11581 50.1 
14376 62.5 14329 5.0 
16858 67.9 16777 56.8 
21500 lager 4 21337 68.2 
27000 91.4 26830 83.1 
32934 100.3 32830 86.9 
Yaw Angle = 80 degrees : 
2849 22.0 2845 rer 4: 
3543 24.2 3553 Day, 1 
5326 30.4 5302 24.6 
7461 39.0 7466 ie Wi | 
11566 49.3 11643 38.8 
14388 56.0 14364 43.0 
16934 61.7 16898 49.4 
21587 69.0 21548 55%. | 
273547 83.7 27427 65.0 
33110 95.3 33075 74.5 


Yaw Angle = 90 degrees 
(Parallel case) 


2849 23.6 2848 7.0 
3548 Pes: 3590 19.4 
5926 S16 a, 5287 Zor) 
7466 43.5 7461 op es) 
11569 Sects) 11609 43.0 
14364 64.6 14310 49.3 
16934 71.8 17018 55.8 
21838 81.9 21591 6225 
27222 oy ag 27180 75.9 
33120 110.4 33214 67.2 
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Table IV. 


Reynolds Number 


1435 
2140 
2727 
3703 


4991 . 


6722 
9001 
11203 
16203 


1435 
2140 
2727 
3703 
4984 
6722 
9001 
Liaz9 
16441 


1435 
2124 
ALT fad 
3703 
4978 
6722 
9019 
Liars 
16444 


Inch Diameter Model 


Basic Model 


Yaw Angle = O degrees 
(Normal case) 


Average Nusselt Number 


Yaw Angle = 10 degrees 


Yaw Angle = 20 degrees 
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Summary of Heat Transfer Data for the 0.25 


197 
2350 
250 
Se 
36.6 
44.1 
50 20 
56.9 
74.9 


20.1 
UBS, 
ZO 
32.4 
Ses 
44.2 
D2 al 
59.0 
76.9 


20.4 
24.6 
ITE, 
S558, 
see 
46.9 
D0. | 
Orel 
79.8 


Table IV. (cont) 
Yaw Angle = 30 degrees 


Basic Model Model with extension 
Re Nu Re Nu 
1435 ZO .2 1415 19.8 
2124 26.5 2145 23.9 
2727 26.9 2684 26.2 
B73 34.7 3740 34.0 
4978 40.6 5028 40.5 
6697 47.9 6751 46.5 
9019 56.2 8975 56.4 

11220 oy5) al 11058 60.8 
16394 80.4 16465 TiO 
Yaw Angle = 40 degrees 
1435 19.6 1415 19.0 
2124 Se | 2145 23.4 
2727 26.6 2684 25.4 
B7O3 33.5 3740 33.0 
4978 40.0 5028 39.4 
6697 47.2 6751 46.5 
9023 56.3 8993 53.4 
12027 62.2 11029 60.8 
16514 (gsveae4 16605 76.8 
Yaw Angle = 50 degrees 
1435 18.8 1415 i7 .8 
2124 22.7 2145 2H ab 
2727 eo 2684 23.8 
3703 SAe,.6 3740 29.9 
4978 37.3 5028 36 .0 
6697 44.0 6751 41.5 
9096 223 8975 51.6 
11187 57.9 10908 Sow 
16689 73.6 16605 7s 
Yaw Angle = 60 degrees 
1435 WieZ 1414 16.0 
2124 2020 2129 122 
2727 22.8 2659 20.8 
3703 293 3740 26", 7 
4991 By 5028 3850 
GHi22 41.6 6775 38.6 
9111 48.5 8975 46.3 
11323 Seer 10969 47.5 
16729 69.2 16720 62.0 
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Taple IV. (Cont) 
Yaw Angle = 7/70 degrees 
Basic Model Model with extension 


Re Nu Re Nu 
1435 14.7 1414 13.8 
2124 iv .3 2129 Lie 
2727 19.6 2659 18.1 
3712 24.7 3740 23.1 
5004 23) 16 5028 27.4 
6722 34.3 6751 3 ig3 
9165 40.1 9048 37.9 

11360 43.2 11089 41.8 

16827 58.5 16779 55m 

Yaw Angle = 80 degrees 
1435 11.7 1414 10.4 

2124 13.6 2129 123 
2727 1 a7 2659 13.4 
3721 18.6 3740 lie 
5004 22.6 5028 202 
6771 26.1 6800 24.5 
9147 g1.9 9085 29a 

11447 B73 11148 32.9 

16729 43.6 16760 42.4 

Yaw Angle = 90 degrees 
(Parallel case) 
1435 10.6 1414 ei! 
2124 12.7 2129 925 
Zo 27 14.7 2671 10.3 
3721 18.3 3740 14.1 
5011 2. © 5028 Liao 
6746 21.8 6800 20.7 
9165 33.0 9103 28.3 
11476 7.2 11206 3230 
16803 SU a8, 16703 44.0 
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Quantity 


Too 


Slope 


Slope 


Table V. 


Uncertainties for Representative runs of the 
0.50 inch Diameter Model 


RUN: B5 
Value Uncertainty 
529.7 0.1 
0.047 0.001 
2098.85 0.052 
542.0 O48 
-5 -5 
0.4226x10 0.0001x10 
-5 -5 
1.2441x10 0.0002x10 
0.00722 0.00005 
14.6 0.2 
3542 38 
37.9 ¢.3 
RUN: El 
529.9 0.8 
0.485 0.010 
2094.01 0.104 
542.2 Ox 1 
-5 -5 
0.4227x10 0.0001x10 
1.24442%107> 0.0002x10~> 
0.0122 0.0004 
NT fa | OF 
11487 118 
64.2 2) 
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(AX) 10 
x 


07012 


O19 
2.060 
0.005 
0.018 
0.024 


0.016 


1.03 


-2 


Table V. 


Quantity 


Too 


Slope 


RUN: I12 
Value Uncertainty 
529.3 0.1 
2.47 0.030 
2114.82 0.312 
541.7 0.1 
0.4224x107> 0.0001x10~? 
1. 2436x1077 0.0002x10~> 
0.0224 0.0016 
105.1 0.6 
25802 157 
Nest 8.4 
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Ax 


)10 


019 


205 


OFS 


018 


024 


016 


orl 


61 


-61 


14 


Table VI. Uncertainties for Representative Runs of the 


0.25 inch Diameter Model 
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RUN: A3 

Quantity Value Uncertainty ( Ax) 10° 
es 530.8 oe! 0.019 
Pa 0.031 0.001 3.226 
P 2106.03 0.052 0.002 
a 543.1 0.1 0.018 
k 0.4234x10. 0.0001x107> 0.024 
Mm 1.2460x107 0.0002x10~> 0.016 

Slope 0.01562 0.00014 0.90 
U_. 11.8 OR2 1.61 
Rey, 1435 23 1.61 
Nu 2084 0.2 0.90 
_ 531.4 Oak 0.019 
Py 0.375 0.01 2.670 
P 2099.26 0.104 0.005 
. 543.7 0.1 0.018 
k 0.4238x10~> 0.0001x10~> 0.024 
lL 1.2470x107> 0.0002x10~> 0.016 

Slope 0.03056 : 0.00136 1 65 
UW hi .2 0.6 193% 
Rey, 4978 67 1.9% 
Nu 40.0 1.8 4.45 


Table VI. (Cont) 


RUN: 17 

Quantity Value Uncertainty ( Ags 1077 
ae 528.6 Omi 0.019 
P, 320 0.08 1.905 
P 2082.93 0.52 0.025 
ie 540.4 0.1 0.018 
k BRAD lsstD = 0.0001x107> 0.024 
L 1.2413x107° 0.0002x107> 0.016 

Slope 0.05266 0.0024 4.56 
Une, 138.0 1.3 0.95 
Re, 16729 159 0.95 
Nu 69.2 30 4.56 
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APPENDIX A 


Data Reduction 

Pure electrolytic copper was selected as the material of the model 
because its high thermal conductivity results in an interior thermal 
conductance which is large in comparison with the exterior convective 
thermal conductance. The interior temperature of the model may be taken 
as uniform at any instant. This is the case of negligible internal re- 
sistance in which the Biot number, the ratio of the internal thermal re- 
sistance of the body to the external resistance of its surface, is very 
small. 

The Biot number Bi =-—=-—— (A-1) 


For this case: 
k >> h and Bi-+0 


An energy balance during cooling of the model yields: 
h A (€(T-T +) = - mc (A-2) 


that is, the rate of energy leaving the model is equal to the rate of 
change of internal energy of the model. Here iB represents the combined 


surface conductance due to convection and radiation, 


A, = ho an he. (A-3) 


where h. is the convective heat transfer coefficient and h. is the ef- 
fective radiation heat transfer coefficient. Because of the low emis- 
Sivity of the polished copper model, h_ is less than one percent of h 
and can be neglected. (See Appendix B). 

The energy balance neglects axial heat conduction. The two copper 


cylinders located at the sides of the model proper provide thermal poten- 


tials equal to that of the model, hence, only the capacitance of the 
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small insulating separator cooling at a different rate will result in 
an axial conduction. An estimate of this amount is less than one per- 
cent and can also be neglected. 

Confidence in these assumptions is obtained in the excellent agree- 
ment of the present data for the normal (unyawed) position with the re- 
sults of previous investigators. 


Neglecting the conduction and radiation terms, equation (A-1) can 


be expressed as: 


UC) h A 


7-1, ~ me aed 


If the model temperature is T. at time 6 and 0, integration of the 


above equation yields the following expression for the temperature-time 


history of the body: 











ha 
a 26 _ = 
To om mc To - To h A 
= =e or ee, ee = 9 (A-5) 
es aie [= te mc 
To ee 
The slope of plot ln =—— versus 9 will be: 
rT = lS 
h A. 
slope = = (A-6) 
so that the average heat transfer coefficient: 
h = slope x = (A-7) 
S 


The mass of the model was determined in a Mettler type H analytic 
balance. This checks closely with the determination of the mass by 
the product of the volume of the model times the density for pure elec- 


trolytic copper. The diameter of the model was measured before and after 
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assembly with the other cylinders. The mass of the model after machined 
down to the final size and polished was determined from the relation 


” -a,(2)° (A-8) 


The properties for copper and bakelite are given in Table II. 


From the model cooling curve, the times corresponding to each 0.05 


millivolt decrease in thermocouple output were taken and a plot In 


— 


T 
CO 
Teg 
versus 8 was made. Figure 8 shows one of such for a typical run. 


The best line was then drawn through the experimental points corres- 


ponding to a thermocouple output, from 0.8 to 0.4 millivolts. 


Evaluation of the Film Temperature 
As mentioned by Baldwin [4], for continuum flow an increase in body 
temperature causes an increase in the heat transfer coefficient. This 


effect led to the use of the "film temperature T_"' for evaluating prop- 


f 


erties in engineering correlations. The film temperature is defined as: 


CS (A-9) 


where Bis, is the free stream temperature indicated for each run by the 


potentiometer and R. is the average wall temperature of the model or 
Te om 
Ww 


ean fin 


WwW 2 


The initial wall temperature T was taken to be that corresponding 
in 
to a temperature difference between the model and the free stream that 


produces 0.8 millivolt of thermocouple output. Similarly the final wall 


temperature TO was taken for a thermocouple output of 0.4 millivolt. 
fin 
In this way, the slope, film temperature and pressure measurements were 


evaluated at a temperature difference between the model and free stream 


corresponding to about 0.6 millivolt of thermocouple output. 
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Evaluation of the Air Density 


The free stream density p.. was evaluated from the perfect gas 


equation: 


ps ee (A-10) 





where P is the absolute upstream static pressure evaluated from the 

manometric reading from a pressure tap located at station No. 1 of the 

test section and the atmospheric pressure indicated by the barometer. 
The film air density was evaluated similarly but using the film 


temperature defined by equation (A-9). 


Evaluation of the Air Velocity 
The air velocity was computed from the dynamic pressure indicated 
by the manometer connected to the pitot-static tube at the upstream 


station of the test section. Bernoulli's equation for a streamline 


2 
ne) ace a 
4 = = 
P 5 P (A-11) 
from which 
O25 
Pa 
oe =k Dp (A-12) 
CO 


Evaluation of the Reynolds Number 


The Reynolds number is defined as 


R, = Bae (A-13) 

It was evaluated in two ways. First using the free stream density D 
CO 

after McAdams [10]. Second using the film air density io) » after 


Hilpert (cited by McAdams). The difference between these two ways of 


evaluation of the Reynolds number was discussed by Collins and Williams 
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[l6é]. For the present tests this amounts to a difference in Reynolds 
number of about two percent. Figure 10 shows the experimental points 
obtained for the 0.50 inches model at yaw angle of 20 degrees using 
both ways of correlation for the Reynolds number. All the results 
given are for the second form of evaluation, i.e., using film density. 
The viscosity and thermal conductivity of air were evaluated from the 


NBS Bulletin [18], at the film temperature. 


Evaluation of the Average Nusselt Number 
The average Nusselt number is the average heat transfer in non- 
dimensional form 


Nu = oe (A-14) 


where k is the air thermal conductivity. 
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APPENDIX B 


Radiation 
For a gray body with an emissivity € » at a temperature T™ 
radiating to an isothermal enclosure at temperature 1h , the nee 


oo 


radiant heat exchange is: 


a= € GO ACT - 4h) (B-1) 


CO 


In terms of the effective radiation heat transfer coefficient, 


= 


r. 
q. = hla (= es (B-2) 
Thus 
aa = =2 2 
ie A Cre eT rs) (B-3) 


Using T = 95°F = 555°R 
= © = ° 
One — 7 FeOrR 


€ = 0.04 (Emissivity of polished copper taken from Eckert [17}) 


(0.04) (0.1714) (107°) (555 + 530) (555° + 530°) 


om 
i} 


0.044 Btu/ 
5 


hr ft °F 


Btu/ 2 


2 eG sec ft’ °F 


This is well below one percent of the average convection heat transfer 


coefficients obtained from the present tests. 
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APPENDIX C 


Sample Calculations 
The following sample calculations are based on data obtained for a 


typical run for the 0.50 inch model. 


Determination of the mass of the model: 


Original diameter (before assembling) D,, = 0.563 in 
Length L = 0.750 in 
Original mass m= 27.2430 gm 
Final diameter (after finished and polished) D = 0.501 in 
Z 2 
' a D " 0.501, _ 
Final mass, m = m, ¢ mie = 27 285 (sey? 21.573 gm 
_ Oe. 508 | 
= ome 0.04755 lbm 
Heat transfer area: 
- _ FW (0.501) (0.750) _ z 
A. = 7(DL = 1A? = 0.00819 ft 

Recorded Data: 

Yaw angle NN = 0 degrees 

Recorded chart Speed = 0.2 in/sec 

Ambient temperature aoe 10% °F 

Free stream temperature too= 0.830 mv. 

Dynamic pressure Pa = 0.495 in-H,0 

Manometric static pressure p = 0.57 in-H,0 

Barometric pressure Die = 29.//7 in.Hga 


Determination of the Film Temperature: 


t =t ae = 0.830 + 0.671 = 1.501 mv. 
in in 
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7 = t + At = 0.830 + 0.450 = 1.280 mv. 
W CoO We. 
fin fin . 


From a chart millivolts versus degrees Fahrenheit prepared for Copper- 


Constantan thermocouples; 


t = 69.9 °F t = 99,3 °F t = 89.7 °F 
CO WwW, W., 
in fin 


ye? 


ct) 


The average wall temperature: t = (t + 
in fin 


e = (99.3 + 89.7)/2 = 94.5 °F 


The film temperature: te = (tt t)/2 = (69.9 + 94.5)/2 = 82.2 °F 


Determination of the air viscosity and thermal conductivity: 


From interpolation of the NBS tables of air properties [18] at a 


film temperature Ty = 542.2°R: 


5 


Viscosity [1 = 1.2444 x 10 ~- lbm/ft sec 


0.4227 x 107° Btu/sec ft °F 


“ Thermal conductivity k 
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Deternnination of the slope: 











T mv. AT, In Al Q@ sec 
Le ES 
1.00 L2c00 0 0 
O795 b2G53 0.051 4.2 
0. 90 Po 0.105 o.2 
Oms5 ele? 02163 13.4 
0.80 1.250 @.223 12.3 
G25 leas 33 On237 2300 
O70 12426 0.356 29.6 
0.65 1.338 0.430 354 
0.60 1.667 0.511 41.9 
O55 1.818 0.598 49.1 
O50 2.000 0.693 56.9 
0.45 AN Gad 0.798 65.4 
0.40 2.500 0.916 ERT 
On) 22657 1.050 85.6 
0.30 352885 1,204 98. 2 
/aryih 
Figure 8 is the plot ln = versus time for this sample run. 
AT 
From Figure 8, 
slope = oe = 0.01222 oe 
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Determination of the average heat transfer coefficient: 


o> Tie (0.04755) (0.0918) _ 2 
h = slope A = 0.01222 0.00819 = (0.00650 Btu/sec ft’. 
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Determination of the average Nusselt Number: 
= h .00650 .501 
aL h D . (0.00650) (0 ann 


k 12 (0.4227%10~”) 





Determination of the density: 


= P + inst.corr. = 29.77 + (-0.13) = 29.64 in.Hga 


atm bar 
P = - p = 13.6 x 29.64 - 0.57 = 402.534 in-H,0 
atm Z 
P = (402.534) (12) (14.696) _ 2094.03 lb£/£t2 
33.9 

7 2094 .03 _ 3 

Pp = (S335) (ezece) = 0-074] Ibm/te 

bh 2094.03 i 3 

iy = (53.35) (542.2) ~ 0.0724 lbm/ft 


Determination of the free stream velocity: 


0.5 0.5 
U = |i Py (P_ - P) ] 2 (14.696) 12 (32.174) Oz 
CO Bro Ss = 73 Cr | (P4/p) 


18.2961 (p,/ po» 


0.495 eee 
0.0741 





18.2961 ( = 47.29 ft/sec 


Evaluation of the Reynolds Number: 


The Reynolds Number correlated with the free stream density after 


McAdams [2] 


_ (47.29)(0.501) (0.0741 x 107°) 


=: = 11757 
12 (1.2444 x 107) 


Re, = (User Pll 
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The Reynolds Number correlated with the free stream density after 
Hilpert (cited by McAdams) 


_ (47.29) (0.501) (0.0724 x 107”) 


= = 11487 
12 (1.2444 x 10°”) 


fe (UP PT 


fa, 


APPENDIX D 


Uncertainty Analysis 
The manner in which the present tests compare with the ideal case 
of the heat transfer for an infinite long cylinder was discussed in 


Section five. 


Uncertainties in physical properties: 


No uncertainties were considered for the thermal properties of 
copper given by Eckert [17]. The values of viscosity and thermal 
conductivity for air given by the NBS bulletin [18] were considered 
correct. Their uncertainty depends only upon the film temperature at 
which they were evaluated. The air density was evaluated at the up- 
stream station using the perfect gas equation. Its uncertainty depends 
upon the upstream static pressure and temperature. Uncertainties in the 


thermal properties of the air were found to be negligible. 


Uncertainty in the mass and lineal dimensions of the model: 


The mass of the model was determined in an analytic balance so that 
the uncertainty in the determination in mass is negligible. The errors 
in lineal dimensions are considered of the order of 0.1 percent and 


neglected. 


Instrumentation inaccuracies: 
In the determination of the uncertainties reported, the method of 
Kline and McClintock [19] was used assuming that all variables have a 


normal Gaussian distribution. 
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Temperature: 

The free stream temperature was obtained by a calibrated copper- 
constantan thermocouple connected to a Leeds and Northrup portable 
potentiometer with an ice reference. Assuming a correct calibration of 
the thermocouple, the estimated possible error is + one half of the 
smallest division of the potentiometer or + 0.0025 millivolts which 


corresponds to 0.1 degrees Fahrenheit. 


Pressure: 

Since the range of pressure varied requiring different manometers 
the uncertainty in pressure varies for different runs. An analysis for 
three representative runs at different Reynolds number was made to allow 
an estimation of the uncertainties at low, middle and high Reynolds 


. numbers. 


Uncertainty in the Nusselt number: 
oe h D 
The average Nusselt number Nu ——— (D-1) 


Substituting the average heat transfer coefficient from equation (A-7): 





_ me slope 
Nu TL 7 (D-2) 
The uncertainty in Nusselt number can be expressed as: 
ANu _ Aslope (D-3) 





Nu slope 


Because the uncertainties in the air thermal properties and the un- 
certainties in the lineal dimensions were estimated as negligible, the 


average Nusselt number depends mainly upon the uncertainty in the plot 
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In(AT, / QT) versus time. The determination of the slope was a graphi- 
cal procedure whose accuracy depended on the scale used and on the value 
of the slope itself. The uncertainty in the slope is greater for the 
faster runs when the times for cooling were small. A time scale from 0 
to 200 seconds, and from 0 to 150 seconds was‘used for the determina- 
tion of the slope of the 0.5 inches and 0.25 inches model respectively. 
The difference in the time scale used accounts for the small difference 


in the uncertainties in the slopes of the two models. 


Uncertainty in the, free stream velocity: 


The free stream velocity can be expressed as: 


_ OFS 
Uy eae, -@, ip) (D-4) 
and using the perfect gas law: 0.5 
a 
ee omer sian (D-5) 


where P is the upstream absolute pressure and Py the dynamic pressure. 
The uncertainty in the atmospheric pressure was estimated as 0.01 in-Hg 
and neglected. 


The expression of the uncertainty in the velocity can be written 














- 0.5 
, 2 2 Ze OS 
Ooo = 1/2 oad At 60 + | AP (D-6) 
Uno Py Vee P 





The uncertainties in temperature and absolute pressure were found to be 
less than 0.1 percent and neglected. With this consideration the un- 


certainty in the velocity becomes: 


Ale (Uosu= A by /2 Py (D-7) 
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depending mainly upon the uncertainty in the dynamic pressure measured 


with the pitot-static tube. 


Uncertainty in Reynolds Number: 


The Reynolds number 
U BB 
co! 
Rey (D-8) 
a 
with the air density evaluated at the film temperature. 


Introducing the free stream velocity from Equation (D-4),the 


Reynolds Number becomes: 


K 








3 Oud) 
Re = ————(p, PT) (D-9) 
H UT, d oe) 
The uncertainty in Reynolds number can be expressed as: 
AT AP At 
R A Z LReZ diez Je fs 
a SEP? + ae ts y+ PB) + Gq’ 
Re mb 
H 
(D-10) 


With the same considerations as before with respect to the uncertainties 
in viscosity, temperature and static pressure, the uncertainty in Reynolds 
number depends basically upon the uncertainty in the free stream velo- 


Cdatsy.. 


a al / aay (D-11) 


A Re, / Re, oo oo 


Tables V and VI present an estimation of the uncertainties for 


representative runs of the 0.50 and the 0.25 inch model respectively. 
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